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It is expected that nucleons and their parity partners become degenerate when chiral symmetry 
is restored. We investigate this question in the context of the thermal transition from the hadronic 
phase to the quark-gluon plasma, using lattice QCD simulations with Nf = 2-1-1 flavours. We 
observe a clear sign of parity doubling in the quark-gluon plasma. Besides, we find that the nucleon 
ground state is, within the uncertainty, largely independent of the temperature, whereas temperature 
effects are substantial in the negative-parity (A*) channel, already in the confined phase. 


I. INTRODUCTION 

The role of discrete and continuous symmetries played 
a fundamental role in the development of the theory of 
the strong interaction. Quantum Chromodynamics. Chi¬ 
ral symmetry breaking and its restoration remain topical 
subjects, mostly due to the creation of the quark-gluon 
plasma at relativistic heavy-ion collision experiments at 
the Large Hadron Collider (CERN) and the Relativistic 
Heavy Ion Collider (BNL). It is expected that chiral sym¬ 
metry will be restored at high temperature, as seen e.g. 
in nonperturbative studies using lattice QCD simulations 

Si- 

In the past decades chiral symmetry restoration at fi¬ 
nite temperature has been studied in great detail in the 
mesonic sector Q . One reason is that mesonic correlation 
functions are relatively easily accessible on the lattice 
[1,@. Moreover, susceptibilities related by chiral symme¬ 
try, such as in the pion and scalar meson channels, can 
now be computed using chiral lattice fermions Un¬ 
like the mesonic sector, the baryonic sector has hardly 
been investigated at finite temperature (early work on 
screening masses from lattice QCD can be found in Ref. 
0, andjin the presence of a small chemical potential, 
in Ref. Q). Nevertheless, understanding the behaviour 
of nucleons in a hadronic medium or in the quark-gluon 
plasma is relevant for heavy-ion collisions, where proton 
spectra are routinely measured and compared to theoret¬ 
ical predictions. Just as for mesons, possible in-medium 
modification of nucleons and other baryons might affect 
signals observed in those experiments. 

In the baryon sector the combination of chiral symme¬ 
try and parity leads to a prediction readily testable in 
QCD: namely that of parity doubling, i.e., a degeneracy 
between channels related by parity, provided that both 
symmetries are realised (the argument will be briefly re¬ 
viewed below). At zero temperature, where chiral sym¬ 
metry is broken, parity doubling is not observed, except 
perhaps in the case of excited hadrons [^. However, 
since chiral symmetry is restored at high temperature, 
it should become relevant in the quark-gluon plasma. 

Recently the question of parity doubling has been 
taken up in Ref. [lOj, where it was studied at three tem¬ 
peratures in quenched lattice QCD. In this paper we 


present what is, as far as we know, the first study of nucle¬ 
ons at hnite temperature in lattice QCD with Nf = 2-1-1 
dynamical quarks, for a range of temperatures below and 
above the deconhnement transition. We find clear indi¬ 
cations of parity doubling, occurring in coincidence with 
the deconhnement crossover. Moreover, within our nu¬ 
merical uncertainty, the mass of the nucleon ground state 
is found to be independent of the temperature of the 
hadronic medium. 


II. NUCLEON PROPAGATION 

The standard interpolation operator for a nucleon, 
which we will consider below, is given by (the material re¬ 
viewed here is well-known, see e.g. the textbooks 0113 ) 

OAr(x,T) = eahcWa(x,T) [u^(x, T)C 75 dc(x, t)] , (1) 

where u, d are the quark helds, a, b, c are colour indices, 
other indices are suppressed and C denotes the charge 
conjugation matrix. Under parity one hnds that 

POAr(x, rjT* = 740 Ar(-x, r), (2) 

and hence operators for the positive and negative parity 
channels are obtained as 

OAr±(x,r) = P±Oiv(x,r), P± = i(l±74). (3) 

We consider the usual euclidean correlators, summed over 
the Dirac indices and projected to zero momentum, 

G'±(t) = J (OAr^(x,T)Ow±(0, 0)) . (4) 

It follows from the properties under euclidean time re- 
hection that, in the case of G+(t), forward (backward) 
propagation in time corresponds to the positive-parity 
(negative-parity) channel. On a lattice at a nonzero tem¬ 
perature T, with 0 < T < 1/T, the negative-parity chan¬ 
nel is then propagating with r_ = 1/T — r. Hence both 
parity channels can be obtained from the same correla¬ 
tor, either G+(r) or G_(t). In the case that the signal 
is dominated by the ground states in both channels, this 
leads to the simple exponential Ansatz, 

G±(t) = A±e-”^±^ -b 


( 5 ) 
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TABLE I. Details of the ensembles. The lattice size is 
Ng X Nt, with the temperature T = l/^arNr). Nc,ig (A^src) 
denotes the number of configurations (sources) used at each 
volume. The spatial lattice spacing is as = 0.1227(8) fm, with 
anisotropy asja-r = 3.5. 


with two masses m±. We note that the spectral repre¬ 
sentation is slightly more complicated than for mesonic 
correlators, and reads 

( 6 ) 

where p±{uj) are the spectral functions in the positive 
and negative parity channels [l3l| . 

Provided that chiral symmetry is unbroken, perform¬ 
ing a chiral rotation on the quark fields immediately leads 
to the result that the two parity channels are degenerate, 
and El 

G±(r) = G^(r) = G±(l/T-r), (7) 

up to overall minus signs. In nature, chiral symmetry is 
broken at zero temperature. Indeed, the ground states 
in the two parity channels differ substantially in mass, 
with the nucleon iV(939) considerably lighter than the 
negative-parity partner iV*(1535). Here we study what 
happens when the temperature is increased and chiral 
symmetry is eventually no longer spontaneously broken. 
We only consider G+(t) and drop the -|- from now on. 


III. LATTICE DETAILS 

We study this question using lattice QCD simula¬ 
tions with Nf = 2-1-1 flavours on anisotropic lat¬ 
tices, with a smaller temporal lattice spacing a,- < a^, 
namely ag/a-r = 3.5. We use a Symanzik-improved 
anisotropic gauge action with tree-level tadpole coeffi¬ 
cients and a tadpole-improved Wilson-clover fermion ac¬ 
tion with stout-smeared links; details on the lattice action 
and parameters at zero temperature can be found in the 
work of the Hadron Spectrum Collaboration El- The 
strange quark mass is tuned to the physical value, while 
the light quarks correspond to Mt^ = 384(4) MeV, with 



FIG. 1. Euclidean correlator G(r)/G(0) as a function of tT. 
At each temperature, the filled symbols indicate the minimum 
of the correlator, within the error. 


M„/Mp = 0.466(3) El- We consider a number of ensem¬ 
bles at nonzero temperature, as summarised in Table HI 
The pseudo-critical temperature Tc is determined via the 
renormalised Polyakov loop. The larger value of Tc than 
expected is primarily due to the light quarks being heav¬ 
ier than in nature. More details on the finite-temperature 
ensembles are available in Refs. El^El- 

For the nucleon operators we employ the form ([3]) , with 
Gaussian smearing for the sources rj and sinks, translated 
over the lattice, using El 

p'= C {1 + nH f 7j, (8) 

where H is the spatial hopping part of the Dirac operator 
and G an appropriate normalisation. Most of the results 
are obtained using k = 8.7 and n = 140 Eli applying 
the same smearing parameters at all temperatures. The 
smearing procedure was tuned to maximise the length of 
the effective-mass plateau in the positive-parity nucleon 
channel at the lowest temperature. The links are APE 
smeared, using one iteration with a = 1.33 El. Correla¬ 
tion functions are generated using the Chroma software 
suite El- ^ discussion of smearing dependence is given 
at the end. More details will be available in Ref. IlSll. 


IV. RESULTS 

The nucleon-nucleon correlators are shown in Fig.[T]on 
a logarithmic scale as a function of rT, for the various 
temperatures we consider (the lowest temperature is not 
shown). In the hadronic phase we observe approximate 
exponential decay for both the forward and the back¬ 
ward propagating mode, but with a clear absence of re¬ 
flection symmetry around tT = 1/2, indicating that the 
positive-parity nucleon, with mass m_|_, is considerably 
lighter than the negative-parity nucleon, with mass to_ . 
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FIG. 2. Euclidean correlator G(r)/G(0) versus t/ut for 
the positive-parity channel (left) and negative-parity channel, 
starting at the opposite side of the lattice (right). 


As the temperature is increased, the correlator becomes 
more and more symmetric; we have used filled symbols 
to indicate at which time slices the correlator is mini¬ 
mal, within the statistical error. This moves towards the 
centre of the lattice, indicating that forward and back¬ 
ward propagation become degenerate, i.e. we find parity 
doubling. 

To find the possible temperature dependence of the 
ground states, we show in Fig. [5] both parity chan¬ 
nels from the same correlator, but now as a function 
of t/ot-, starting at the opposite side of the lattice for 
the negative-parity channel. The lowest temperature is 
now also included. Note that we show the correlator up 
to its minimum. On the positive-parity side, we observe 
considerably less temperature dependence than on the 
negative-parity side. This is further demonstrated in Ta¬ 
ble II, where we show the results from fits to the exponen¬ 
tial Ansatz below Tc- The error combines estimates 
of statistical and systematic uncertainties in quadrature. 
At the lowest temperature, we note that both masses are 
higher than in nature (see also Refs. [H, [l^l), presum¬ 
ably due to the quarks being too heavy, but that the mass 
splitting TO- — TO+ ~ 700 MeV is of the right order, albeit 
with a large error. The positive-parity mass is, within the 
error, temperature independent, while to_ shows signif¬ 
icant temperature dependence. We found that above Tc 
simple exponential fits are no longer reliable. This may 
be due to the absence of a clear ground state, i.e. the 
nucleon is no longer a well-defined particle. 

To study the onset of parity doubling in the confined 
phase, the final column contains the dimensionless ratio 


A = 


TO- — TO+ 
TO- + TO+ ’ 


(9) 


with estimates of statistical and systematic errors added 
in brackets. Note that in nature, A = 0.241 at zero 


T/Tc artrijc arm- 

m+ [GeV] 

m- [GeV] 

A 

0.24 0.213(5) 0.33(5) 

1.20(3) 

1.9(3) 

0.209(28)(082) 

0.76 0.209(16) 0.28(3) 

1.18(9) 

1.6(2) 

0.138(29)(130) 

0.84 0.192(17) 0.28(2) 

1.08(9) 

1.6(1) 

0.197(39)(054) 

0.95 0.198(25) 0.22(4) 

1.12(14) 

1.3(2) 

0.052(35)(190) 


TABLE II. Results from exponential fits to Eq. below Tc, 
with statistical and systematic errors added in quadrature. 
The final column contains A = (m- — m+)/(m- -|- m+) with 
estimates of statistical and systematic errors added in brack¬ 
ets. In nature, A = 0.241 at T = 0. 



FIG. 3. Ratio R{t) in Eq. (1101) as a function of tT. 


temperature, while in the parity-degenerate case, A = 
0. We observe a reduction of A as the temperature is 
increased, with a rather large systematic error arising 
mostly from the difficulty of determining the groundstate 
mass in the negative-parity channel at finite temperature. 

In order to investigate parity doubling in more detail 
also above Tc, we consider the following ratio of correla¬ 
tion functions 0 


R{t) 


G{r)-G{l/T-T) 

G(T)+G(l/T-r)- 


( 10 ) 


Note that R{l/T — t) = —R{t) and hence R{1/2T) = 0. 
Consider first low temperature. In the (extreme) case 
that the correlator is dominated by the positive-parity 
ground state and that the negative-parity ground state 
is much heavier, to- ^ m+, we find that R{t) = 1, 
except near tT =1/2, due to the reflection asymmetry. 
On the other hand, in the case of parity doubling, with 
G(t) = G(l/T — r), we find that R{t) = 0. Hence this 
ratio lies naturally between 0 and 1. 

The results for i?(r) are shown in Fig. [31 We ob¬ 
serve that the ratio is distinctly different from zero in 
the hadronic phase, indicating the absence of parity 
doubling, with an appreciably lighter ground state in 
the positive-parity channel. The ratio is fairly constant 
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FIG. 4. Averaged ratio R in Eq. m as a function oi T/Tc- 
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FIG. 5. As in Fig. [4] using point sources and different levels 
of smearing. 


nonzero value below Tc to a value close to zero above Tc- 
The transition coincides surprisingly well with the tran¬ 
sition to the deconfined phase, which is based on the be¬ 
haviour of the Polyakov loop [i3- Since deconfinement 
and chiral symmetry restoration are expected to occur 
around the same temperature Qi, the observed parity 
doubling can hence be explained through the restoration 
of chiral symmetry in the quark-gluon plasma. 

In order to achieve these results, we found that it is 
essential to suppress excited states at early Euclidean 
times, either by using smeared sources and sinks, or by 
considering only a restricted time interval. This is in par¬ 
ticular pertinent for the Wilson-clover lattice fermions 
we use: since the Wilson term violates chiral symmetry 
and becomes relevant at larger energy scales, parity dou¬ 
bling of higher excited states cannot be expected. This is 
demonstrated in Fig.jSJ where we show the dependence of 
R on various sources and sinks. Point sources (no smear¬ 
ing) couple strongly to short-distance states, for which 
chiral symmetry remains explicitly broken by the Wil¬ 
son mass term. This can be partially taken into account 
by excluding the first few time slices in the analysis and 
considering only t/o,- > 5. Applying a small number of 
smearing steps (n = 10), but preserving all time slices, 
has approximately the same effect, while also excluding 
the early-time interval improves the signal considerably. 
Finally, applying more smearing steps (n = 60,140) al¬ 
lows for the signal of parity doubling to fully emerge. 
We note that since chiral symmetry breaking for Wil¬ 
son fermions at large energy scales is a lattice artefact, 
it would be interesting to repeat this analysis using chi- 
rally symmetric fermions (see e.g. Ref. [2^ for a study of 
parity partners using domain-wall fermions at zero tem¬ 
perature), especially in the chiral limit. Finally, we have 
also verified that using different sources (with the same 
quantum numbers) yields qualitatively the same result 


across the entire euclidean time range: the drop towards 
zero as tT approaches 1/2 follows from the symmetry of 
i?(r), as explained above, while for small tT the effects 
of excited states enter. As the temperature is increased, 
the ratio decreases nearly uniformly across the entire time 
range and is close to zero in the quark-gluon plasma, the 
signal of parity doubling. 

In order to quantify the parity degeneracy further, we 
consider the average ratio i?, defined as 

^ -N —1 ’ ^ 

where (t(t) denotes the error in R(t) and t„ = uQt- 
Again, the natural value of R is close to but below unity 
in the chirally broken phase and close to zero in the 
parity-doubling phase. The results for R are shown in 
Fig. m We observe a clear crossover behaviour from a 


V. SUMMARY 


We have carried out a study of nucleons at finite tem¬ 
perature using lattice QCD simulations with Nf = 2 -|- 1 
dynamical flavours over a range of temperatures. We 
found approximate parity doubling, intimately linked 
with the transition to the deconfined phase. We found 
that the positive-parity nucleon mass is largely indepen¬ 
dent of the temperature in the hadronic phase, while the 
negative-parity channel shows clear temperature depen¬ 
dence already below Tc- In the context of heavy-ion colli¬ 
sions, we note that the temperature independence of the 
nucleon justifies treating the proton as unmodified by 
the thermal medium, while a temperature-dependent N* 
mass may have implications for heavy-ion phenomenol¬ 
ogy [Hi. 

As an outlook [T3l| . we note here that we are currently 
considering baryons containing strange quarks, with the 
aim of determining the effect from the explicit breaking 
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by the larger strange quark mass. We note that the role 
of excited states will be of particular interest when com¬ 
paring nonchiral fermions, such as in this paper, with 
chiral fermions, such as domain wall fermions. Finally, it 
will also be interesting to analyse our results in terms of 
baryon spectral functions, using Eq. ([6]), and to address 
them in terms of effective models, such as the one pro¬ 
posed in Ref. 2M (see also Ref. for a lattice study of 
the axial charge of negative-parity nucleons in this con¬ 
text). 
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